#include "functions.h"

int main()
{
  for(int k = 5; k < 21; k++)
    {
      double** H = new double*[k];
      for(int i = 0; i < k; i++)
	{
	  H[i] = new double[k]();
	  for(int j = 0; j < k; j++)
	    H[i][j] = 1.0/(i + j + 1.0);
	}
      std::cout << "The Condition Number of " << k
		<< "-order Hilbert Matrix is: ";
      std::cout << Condition_Number_Norm_Infinite(H, k) << std::endl;
      for(int i = 0; i < k; i++)
	delete [] H[i];
      delete [] H;
    }
  return 0;
};
